Scaling relations in the collapse of elastic icosadeltahedral shells under uniform 

external pressure 
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We discuss collapse of icosadeltahedral shells subjected to uniform external pressure. We demon- 
strate that there is a universal collapse pressure curve. The parameter that uniquely determines 
the collapse pressure is shown to be the Foppl-von Karman number of the non-pressurized shells. 
Numerical results are interpreted in terms of buckling instabilities of spheres and cylinders under 
hydrostatic pressure. 
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It has been experimentally demonstrated that micro- 
scopic and nanoscopic capsules can collapse under ex- 
ternal pressure, evaporation of solvent, or point-forcing 1 . 
In most of these cases, the studied capsules were nearly 
spherical. Viruses represent an exception as they can 
take on different shapes. Even icosahedral viruses that 
are studied in this article can have very different shapes, 
some of them being nearly perfectly spherical, while oth- 
ers have pronounced polyhedral shape and nearly flat 
faces^. These viruses can be represented as triangulated 
shapes of spherical topology with an icosa(delta)hedral 
backbone - this is known as the principle of quasiequiv- 
lence and forms the basis of the Caspar-Klug classifica- 
tion of icosahedral viruses 3 -. This also means that viruses 
can be thought of as consisting of clusters of proteins, 
with 10(T — 1) of these clusters consisting of six pro- 
teins (hcxamers) and exactly twelve of them consisting 
of five proteins (pentamers). Here T is the triangulation 
number of the virus and describes the arrangement and 
number of protein subunits on the capsid 4 * 5 -. 

The central issue of this work is to explore whether and 
how the equilibrium shape of an icosadeltahedral shell 
(viral capsid in particular) influences its response to the 
external mechanical pressure. This should be important 
for empty viral capsids submitted to osmotic pressure 
of the external solution. The impenetrability of the vi- 
ral capsid to the osmoticant, such as polyethylene glycol 
(PEG), gives rise to a mechanical pressure across the 
capsid shell that compresses it. Experiments along these 
lines on complete virions (capsids as well as genome) are 
performed on bacteriophages in PEG bathing solutions 11 . 
Our interest here is not in the pressure dependence of the 
amount of DNA encapsidated in these viruses 1 ^, but in 
the response of empty capsids under mechanical pressure 
and, most importantly, their collapse. 

By applying the continuum elasticity theory to viral 
capsids (empty viruses), Lidmar et a.' 4 (LMN in the fol- 
lowing) have shown that their shape can be understood 
in terms of a single parameter, the so-called Foppl-von 
Karman (FvK) number (7) given as 



j = Y(R) 2 /k. 



(1) 



Here Y, and k are the two-dimensional Young's modulus 
and the bending rigidity of the viral protein sheet, respec- 
tively, and (R) is the mean radius of the viral capsid. For 
7 smaller than about 250, the equilibrium shape of the 
capsid that minimizes its elastic energy is practically a 
perfect sphere [see Fig. [lb)]. When 250 < 7 < 5000, 
a continuous transition in equilibrium shape takes place 
and the capsids assume a more aspherical shape. LMN 
have termed this the "buckling transition" since the re- 
gions surrounding the pentagonal disclinations (protein 
pentamers) "buckle out" from the sphere so that the sur- 
face surrounding each of the disclinations is nearly coni- 
cal (it should be understood that the buckling transition 
in the LMN terminology is not a consequence of any ex- 
ternal forcing of the capsid but simply results from the 
minimization of elastic energy of the shell) . In the region 
250 < 7 < 10 4 , the capsids can be imagined as a union of 
twelve conical frusta (with apices at the icosahedron ver- 
tices) that are fastened together at their bases 5,6 . This 
approximation provides a good account of the energet- 
ics of the shells 4,5,7 . Within the range 10 4 < 7 < 10 6 , 
the conical description of the shell becomes less satis- 
factory with regards to the shell energetics [c.f. Fig. 6 
of Ref. i and Fig. []}>)] since another creeping transi- 
tion takes place that flattens the icosahedron faces and 
at the same time sharpening the regions around their 
edges. This effect has been explained by Witten and Li£ 
and Lobkovsky^ as originating from the stretching en- 
ergy along the edges that becomes prohibitively large as 
the shell size increases (or as 7 increases) 7 -. The scaling 
relations that are characteristic for the ridge sharpen- 
ing are observed^ when 7 > 10 6 , so that in the region 
10 4 < 7 < 10 6 neither cones nor ridges can provide ad- 
equate representation of the shell shape and energetics. 
Note, however, that the borders of different regimes of 
the FvK number are quite smeared, especially towards 
the ridge sharpening regime, since the transitions in shell 
shapes are continuous. As the FvK number scales with 
(i?) 2 , larger viruses should look more polyhedral, and 
this is indeed a general trend observed in the experimen- 
tal data (note, however, that the influence of the viral 
genome on the shell shape has been neglected, so that, 
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strictly speaking, the above consideration should be ap- 
plied only to empty viral capsids 1 ^). 

Since the equilibrium shape of the capsid/shell changes 
with the FvK number, we could a priori expect that the 
response of the shells to external pressure is also strongly 
influenced by its FvK number. An intriguing question to 
pose is whether the FvK number uniquely determines 
the response of the shells to external pressure. This has 
in fact been proven in Ref. [TH where the authors con- 
sidered the susceptibility of shells as a function of their 
FvK number. That study has, however, concentrated on 
elastic response of the shells to external forcing whereas 
we shall be interested in the collapse of shells. Interest- 
ingly, even non-pressurized empty viral shells have been 
found to be prone to collapse in two recent molecular 
dynamics studie a 14 ' 15 . Our study should be in many re- 
spects complementary as it allows a study of the collapse 
as a function of elastic (coarse grained) parametrization 
of the protein-protein interactions in viral capsids. 

We model the shell as a polyhedron with the icosahe- 
dral order whose neighboring vertices are connected with 
springs, so that the stretching energy is given by 

where e is the spring constant, a is the equilibrium sepa- 
ration of the neighboring (i and j) vertices, and is the 
vector pointing at the i-th vertex. The bending energy 
of the shell is given by 

^ = |E( nj - nj ) 2 , (3) 

i,j 

where / and J are two polyhedron faces sharing an edge 
and nj (nj) is the unit normal of the 7-th (J-th) face. 
The macroscopic elasticity constants of the shell material 
can be derived from e and k as demonstrated in Ref. |4j, 
so that the Young's modulus is Y = 2e/\/3, the (mean) 
bending rigidity k — \^3k/2, the Gaussian bending rigid- 
ity kq = — 4«/3, and the Poisson ratio v = 1/3. The 
shape of the shell under external pressure (p) is found by 
minimizing the Hamiltonian of the problem, 

H = H s + H b +pV, (4) 

where V = J2i r -f,i ' ( r J,2 x r /,3)/6 is the shell volume, and 
r i,i) r 7,2i and are the vertices of the I-th triangular 
face in the clockwise order. For given shell T-number and 
parameters a, e, and k, the non-pressurized shape of the 
shell is uniquely determined by its FvK number. Upon 
increasing pressure, the minimum-energy shape of the 
shell [minimizing H in Eq.|4])] changes, i.e. the shell de- 
forms. We perturb the vertices of the thus obtained shape 
by adding a random displacement to each of the vertices, 
v' i = Vi + uei, where ej is a random three-dimensional 
unit vector, different for each vertex, and u is the am- 
plitude of the displacement. The perturbed shape can 
be thought of as a particular conformation excited by 



(low-, h < a) temperature fluctuations and the pertur- 
bation has in general non-vanishing projections on each 
of the vibrational eigenmodes of the shell (see below). 
We then again minimize the energy of the perturbed 
shape. For sub-critical pressures, the shape equilibrates 
back to the unperturbed state, but at some critical pres- 
sure (p c )j the new conformation that is adopted by the 
shell differs from the unperturbed state, i.e a discontinu- 
ous transition in the shell shape parameters and volume 
is observed. Note here that we assume that the elastic in- 
stability of the shell takes place prior to its rupture which 
is corroborated by the experiments on empty bacterio- 
phage capsids which are shown to withstand large forces 
and indentations by the tip of an atomic force micro- 
scope (~ 6 nm) prior to nonlinear response and possible 
ruptures*. To relate our results to those previously pub- 
lished for non-pressurized icosadeltahedral shells^^ and 
continuum theories of shell instabilit y 17 ' 18 ! 19 , we shall be 
particularly interested in shells with large T-numbers. 
The results should be therefore directly applicable to 
viruses with large number of protein subunits, although 
the numerical simulations can be even more easily per- 
formed for shells with small T-numbers. 

Deformations of elastic plates subjected to external 
forces are known to belong to two distinct categories, 
depending on their elastic parameters. It has been exper- 
imentally demonstrated that deformations of a clamped 
half-cylindrical surface under point forcing 20 and the 
gravity induced draping of naturally flat, isotropic 
sheet^i can be understood in terms of the creation of con- 
ical singularities (the so-called d-concs) 22 . This means 
that the energy of deformation is dominantly of the bend- 
ing type. On the other hand, the crumpling of a sheet of 
paper has been described as dominated by the generation 
of narrow stretching ridge o 23 ' 24 . In this case, the ener- 
getics of deformation is still dominated by the bending 
contribution, but the stretching energy becomes of com- 
parable magnitude 9 . As already discussed, in the shapes 
and energetics of equilibrium non-pressurized shells, both 
types of deformations can be seen but in different regimes 
of the FvK number. Thus, we expect that the same two 
types of deformations should also be seen in the energet- 
ics of the shell collapse and this is one of the central issues 
of this work. LMN have demonstrated^ that the elastic 
energy of non-pressurized shells scales as E cx Kf(j), 
where the function 7(7) assumes different forms depend- 
ing upon whether the energy of the shell is dominated by 
stretching or bending (i.e. on the value of 7). A simple 
analysis suggests that if there is universality in the col- 
lapsing pressure, and if it is related to non-pressurized 
energetics and shape of the shell, it should reveal itself in 
the scaling p' c — ► p c (R) 3 /n, where (R) is the mean radius 
of the non- pressurized shell. The results presented in Fig. 
0J,) convincingly demonstrate the universality in p' c . This 
figure displays the scaled critical pressures p' c as a func- 
tion of FvK number of the shell in the non-pressurized 
state. Note that the thus rescaled collapsing pressures 
for shells of different elastic properties (Y and k) and 
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T-numbers (i.e. (R) if a is fixed) all fall on the same 
universal curve which we denote by U (7) . The collapsing 
pressure is thus given as 



Pc 



(R) 



Mi), 



(5) 



Intriguingly enough, for values of 7 that should be typ- 
ical of large viruses (300 < 7 < 3000) 4 -, U{^) is nearly 
constant, i.e. U(^f) m 70. The measurement of collapsing 
pressures should thus produce a direct information on k. 




FIG. 1: Panel (a): Scaled critical pressures (symbols) of the 
icosadeltahedral shells with e = 1 and T = 441 (+), e = 10 
and T = 441 (x), e = 5 and T — 100 (squares), e = 5 
and T — 49 (circles), and e = 4 and T = 73 (triangles) 
as a function of the FvK number 7. In these calculations, 
u = 0.005a, e was kept fixed and « so to produce a variation 
in 7. The dashed lines show scalings with 7 as discussed 
in the text. Panel (b): The energetics of non-pressurized 
shells as a function of 7. The regimes in which the shell can 
be represented as a sphere, assembly of cones or ridges are 
indicated. Analytical expressions for the shell energetics are 
indicated by full (sphere) , dashed (cones) and dotted (ridges) 
lines^. 

The collapse of shells as the pressure increases can be 
understood by following the eigenmodes of shells for a 
given pressure. These are obtained by diagonalization 
of the matrix of second derivatives of the Hamiltonian 



in Eq. |QJ, viz. d 2 H/dridrj. Near the collapsing pres- 
sure, the frequency of one of the eigenmodes approaches 
zero indicating the instability of the shell with respect 
to the motion pattern specified by the eigenvector of the 
critical eigenmode. This approach has been successfully 
applied in the analysis of ridge buckling instability- 5 . 
We have found that the eigenvector patterns of critical 
eigenmodes change with 7, so that different shell mo- 
tions lead to its instability, depending on 7. One could 
proceed further with this analysis but a lot of insight 
in the numerically obtained results can be obtained by 
applying existing analytical results for the critical hydro- 
static pressures leading to buckling of sphere o 17 ' 18 ! 19 and 
cylinders 18 . Namely, although the shapes of both pres- 
surized and unpressurized shells are fairly complex, their 
surface can be separated in regions that locally resemble 
spheres (two finite radii of curvature) or cylinders (one 
radius of curvature tending to infinity) , depending on the 
FvK number. We thus expect that the scaling of criti- 
cal pressures for spheres and cylinders should be seen 
in our numerical data in different regions of 7. For FvK 
number smaller than about 250, the equilibrium shape of 
the unpressurized shell is practically a perfect spher o 4 ' 10 . 
Therefore, the critical pressures should be proportional 
to 



17.18.19 



^sphere K 



"Y7 



(6) 



which is an expression for collapse pressure of thin 
spheres. This is exactly what we find from our numerical 
results in the region 90 < 7 < 250 as can be seen by 
comparison of a function proportional to 7 1 / 2 with the 
data. A deviation from this behavior occurs for 7 < 90, 
but note that for a shell made of uniform elastic material 
7 = 12(1 — v 2 )(R/d) 2 , where d is an effective thickness of 
the shell. Thus, for 7 < 90, R/d is no longer negligible 
and an expression for thin shells in Eq. ^ becomes in- 
applicable. For shells that underwent the buckling tran- 
sition in their non-pressurized state, the shell surface is 
nearly conical around the pentagonal disclinations (pro- 
tein pentamers) i.e. locally cylindrical, so that the critical 
pressure should be proportional to^ 8 - 



cylinder 
Pc (x 



R 3 : 

U cyl 



(7) 



where R cy i is the maximum radius of the cylindrical re- 
gions in the shell (this suggest that the critical points of 
collapse are located in the middle of icosahedron edges). 
The maximum radii of the cones are R cy i ~ (R) , so that 
the critical pressures should be independent of 7, which is 
also obtained in our numerical results for 300 < 7 < 3000 
(i.e. U("j) cx const.). As the FvK number increases, the 
shell shape just prior to collapsing shows flattening of the 
faces and concentration of curvatures along the ridges, so 
that the radius of cylindrical regions immediately before 
collapse notably decreases from its non-pressurized value 
as a result of applied pressure. This can also be seen 
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as an increase in the collapse pressure, i.e. its deviation 
from the constant value predicted by Eq. ([7]). Note also 
that the description of shells in terms of cones becomes 
inaccurate when 10 4 < 7 <■ 10 6 [see Fig. [1J>)] and nei- 
ther ridges nor cones provide an adequate description of 
the shell shape in this region of FvK number since the 
curvature along ridges scales in a complicated way with 
R and 7. When 7 > 10 6 , the ridge-sharpening transition 
gradually takes place, and the characteristic curvature 
radii of (non-pressurized) ridges scale as^ 



R, 



cyl 



(R)r 



(8) 



which means that the critical pressures for shell collapse 
should scale as 



Pc oc 



i? 3 



«V7 

(JZ)3' 



(9) 



i.e. functionally the same as for shells that are spheri- 
cal in their non-pressurized state [7 < 250, see Eq. ©]. 
The scaling of p c with 7 1 / 2 in this regime (^(7) oc 7 1 / 2 ) 
is again strongly suggested by our numerical results. A 
detailed exploration of this region of large FvK number 
requires numerical studies of very large shells (large T- 
numbers) since the lower bound on radii of ridge cur- 
vatures is set by a. This slows down the calculations, 
but from our results we estimate that T = 441 shells 
should be a reliable model for continuum shells to about 
7 ~ 10 8 . The region of such large 7 is not important for 
viruses but may be of importance for other shell struc- 
tures with icosahedral order such as e.g. self-assembled 
hollow icosahedra in salt-free catanionic solutions 26 . 



We have confirmed that the main features of our nu- 
merical results on the collapse of elastic icosadeltahedral 
shells under uniform external pressures can be under- 
stood on the basis of buckling instability pressures that 
are characteristic of spherical and cylindrical surfaces. 
Note that the scaling obtained in our results is also pre- 
dicted by Eqs. ©, 0, and ©. The knowledge of bend- 
ing rigidity of viral capsids is required in order to apply 
our results to viruses. This quantity is unfortunately 
not known precisely, but can be estimated to be of the 
order of k ~ 40/cbT from a previous study of bacterio- 
phage shells under internal pressure generated by DNA 10 
(a lower value of about 10 - 15 fcsT has been found in 
Ref. H3) ■ Assuming a radius of about (R) ~ 30 nm that 
is typical for e.g. A-bacteriophage produces collapsing 
pressures in the regime of FvK numbers that are typ- 
ical for large viruses of about p c ~ 5 atm. The pres- 
sures of this order of magnitude can be easily achieved in 
experiments 1 -^, but we again stress that a test of our pre- 
diction would require exposing empty viruses to osmotic 
pressure. In our elastic model we have assumed that the 
spontaneous curvature of the viral shell is zero. While it 
is still a matter of debate whether the protein coatings of 
viruses do indeed have some non-vanishing spontaneous 
curvature, this is suggested by a recent study^ 8 .. This 
may also have some influence on the precise value of col- 
lapse pressures. 
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